Table1401 <- read.csv("http://www.statpower.net/Content/311/Lecture Notes/Table1401.csv")
Table1401
library(reshape)
rm.data <- melt(Table1401,id=c("Subject"),measured=c("Depression"))
colnames(rm.data)[2:3] <- c("Season","Depression")
head(rm.data,20)
library(ez)
results <-ezANOVA(data=rm.data,dv=.(Depression),
    + wid=.(Subject),within=.(Season),type="III",return_aov=TRUE)
anova(lm(Depression ~ factor(Subject) * Season, data = rm.data))
Table1503 <- read.csv("http://www.statpower.net/Content/311/Lecture Notes/Table1503.csv")
names(Table1503)
Subject<-1:6
Table1503 <- cbind(Subject,Table1503)
temp <- melt(Table1503,id.vars=1,measure.vars=2:10)
A <- rep(c(rep(1,6),rep(2,6),rep(3,6)),3)
B<- c(rep(1,18),rep(2,18),rep(3,18))
rm.data <- data.frame(cbind(temp$Subject,temp$value,A,B))
colnames(rm.data) <- c("Subject","Rating","Orientation",
                        "Distortion")
rm.data$Subject <- factor(rm.data$Subject)
rm.data$Orientation <- factor(rm.data$Orientation)
rm.data$Distortion <- factor(rm.data$Distortion)
head(rm.data,18)
anova(lm(Rating~Subject*Orientation*Distortion,data=rm.data))
ezANOVA(rm.data,wid=.(Subject),dv=.(Rating),within=.(Orientation,Distortion))
#############
Table1505 <- read.csv("http://www.statpower.net/Content/311/Lecture Notes/Table1505.csv")
rm.data <- data.frame(melt(Table1505,id.vars=1:2,measure.vars=3:6))
rm.data$A <- factor(rm.data$A)
colnames(rm.data) = c("Subject","Method","Time","Score")
head(rm.data,18)
ezANOVA(rm.data,wid=.(Subject),dv=.(Score),between=.(Method),within=.(Time))